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Abstract A massive intruder in a homogeneously driven granular fluid, in 
dilute configurations, performs a memory-less Brownian motion with drag and 
| temperature simply related to the average density and temperature of the 

fluid. At volume fraction ~ fO — 50% the intruder's velocity correlates with the 
local fluid velocity field: such situation is approximately described by a system 
of coupled linear Langevin equations equivalent to a generalized Brownian 
motion with memory. Here one may verify the breakdown of the Fluctuation- 
Dissipation relation and the presence of a net entropy flux - from the fluid to 
the intruder - whose fluctuations satisfy the Fluctuation Relation. 
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Granular fluids represent a valid benchmark for modern theories of non- 
q equilibrium statistical mechanics [3]. Due to dissipative interactions among 

the microscopic constituents, energy is not conserved and an external source 
is necessary to maintain a stationary state. The consequence is a breakdown of 
time reversal invariance and the failure of properties such as the Equilibrium 
Fluctuation-Dissipation relation (EFDR) [BJ. In recent years, a systematic the- 
ory for the dilute limit has been developed, in good agreement with numerical 
simulations [2JH] , while a general understanding of dense granular fluids is still 
lacking. A common approach is the so-called Enskog correction [2j|3], which 
reduces the breakdown of Molecular Chaos to a renormalization of the col- 
lision frequency. In cooling regimes, the Enskog theory may describe strong 



o 

(N 



All authors 

CNR-ISC c/o Dip. di Fisica, 
Universita' degli Studi di Roma " La Sapienza" , 
Piazzale Aldo Moro 2, 1-00185 Roma, Italy. 
Tcl.:+39-0649913508 Fax:+39-064463158 



H ' E-mail 



E-mail 
E-mail 
E-mail 



andrea.puglisi@romal.infn.it 
ale.sarracino@gmail.com 
ggradcnigo@gmail.com 
villamaina@gmail.com 



2 



A. Puglisi et al. 



non-equilibrium effects, due to the explicit cooling time-dependence [32]. Nev- 
ertheless it cannot describe dynamical effects in stationary regimes, such as 
multiple characteristic times or different decays of response and autocorrela- 
tion 0HH]. 

Here we review a recent model (24) for the dynamics of a massive tracer 
moving in a gas of smaller granular particles, both coupled to an external 
bath. Taking as reference point the dilute limit, where the system has a closed 
analytical description [23] , a Langevin equation linearly coupled to a fluctuat- 
ing local velocity field is proposed as first approximation capable of describing 
the dense case. Its main features are: i) the decay of correlation and response 
functions is not simply exponential and shows backscattering [15114] and ii) 
the EFDR pHl[T5] of the first and second kind do not hold. In such a model, 
detailed balance is not necessarily satisfied, and a fluctuating entropy produc- 
tion [35] can be measured, which fairly verifies the Fluctuation Relation [TT] 

H2E3]. 

The model reviewed here is the following: an "intruder" disc of mass mo = 
M and radius R, moving in a gas of N granular discs with mass m, = m 
(i > 0) and radius r, in a two dimensional box of area A — L 2 . We denote 
by n — N / A the number density of the gas and by <fi the occupied volume 
fraction, i.e. <j> = Tr(Nr 2 + R 2 )/A and we denote by V (or vo) and v (or Vi with 
i > 0) the velocity vector of the tracer and of the gas particles, respectively. 
Interactions among the particles are hard-core binary instantaneous inelastic 
collisions, such that particle i, after a collision with particle j, comes out 
with a velocity v\ = t>j — (1 + a) m " J m [(t>, — Vj) ■ n]n where h is the unit 
vector joining the particles' centers of mass and a £ [0, 1] is the restitution 
coefficient (a = 1 is the elastic case). The mean free path of the intruder 
is proportional to Iq = l/(n(r + R)) and we denote by t c its mean collision 
time. Two kinetic temperatures can be introduced for the two species: the gas 
temperature T g = m{v 2 )/2 and the tracer temperature T tr = M{V 2 )/2. 

The equation of motion of the i-th particle reads: rriiVi(t) = —jbVi(t) + 
fi(t)+£b{t). Here fi(t) is the force taking into account the collisions of particle 
i with other particles, and £b(i) is a white noise (different for all particles), 
with (&(*)) = and (&(*)&(*')) = 2T b ~/ b 6{t - t'). The effect of the external 
energy source balances the energy lost in the collisions and a stationary state 
is attained with m^vj) < T b ^MMIMMi ■ 

At low packing fractions, <f> < 0.1, and in the large mass limit, m/M <C 1, 
using the Enskog approximation it has been shown [23] that the dynamics of 
the intruder is described by a linear Langevin equation: 

MV = -r E v + £ E , (1) 

with S E a white noise with (£ E ) = 0, {£ E {t)£ E {t')) = 25(t-t')r E Tg andT t f = 
(7&Tf, -(- r )gl^ L Tg)/r E is the tracer's temperature. In this limit the velocity 
autocorrelation function shows a simple exponential decay, with characteristic 
time M/T E , where r E = - lb + 7^ and 7^ = 92( \+ R) ^/2^mT g {\ + a) where 
9i(r + R) is the pair correlation function for a gas particle and the intruder 
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at contact. Time- reversal and the EFDR, weakly modified for uniform dilute 
granular gases [T71I51I21]. become perfectly satisfied for a massive intruder. 

As the packing fraction is increased, the Enskog approximation fails in pre- 
dicting dynamical properties. In particular, velocity autocorrelation C(t) = 
(V(t)V(0))/(V 2 ) and linear response function R(t) = SV(t)/SV(0) show an 
exponential decay modulated by oscillating functions [4"ll24| . Moreover viola- 
tions of the EFDR C(t) = R(t) are observed for a < 1 [151125] . The Enskog 
approximation is unable to explain the observed functional forms, because it 
only modifies by a constant factor the collision frequency [2,23 : a model with 
more than one characteristic time is needed. A first approximation is given by 
an auxiliary field coupled to the intruder's velocity: 

mv = -r E (v -u) + ^JW7r g e v (2) 

M'U = -r'U - r E V + ^/2TT b 8u, 

where By and £ u are white noises of unitary variance. Two new parameters 
appear: the mass of the local field M' and its drag coefficient r' . The dilute 
limit here is obtained for r' ~ M' — > oo. In such a limit indeed U — >■ and 
the equation for V comes back in the form discussed above [33]. In such a 
form $2§, the dynamics of the tracer is remarkably simple: indeed V follows 
a Langevin equation in a Lagrangian frame with respect to a field U, which 
is the local average velocity field of the gas particles colliding with the tracer. 
A first justification of this model comes from realizing [21] that it is equiva- 
lent to a Generalized Langevin Equation with exponential memory, which is 
consistent with a typical approximation done for Brownian Motion when, at 
high densities, the coupling of the intruder with fluid hydrodynamic modes, 
decaying exponentially in time (see [28], Cap. 8.6 and 9.1), must be taken into 
account. Such a coupling, which in principle involves a continuum of modes, is 
reduced here to a single dominant mode: this is sufficient to introduce a new 
non-trivial timescale. The full coupling would reproduce finer features which 
become relevant at larger densities or larger times, such as long-time power- 
law tails. The fact that the "temperature" of the local velocity field U is equal 
to the bath temperature Tf, comes as a consequence of the conservation of 
momentum in collisions, implying that the average velocity of a group of par- 
ticles is not changed by collisions among themselves and is only affected by the 
external bath and a (small) number of collisions with outside particles. This 
scenario is fully consistent with recent study of hydrodynamic fluctuations for 
the velocity field of the same fluid model [ZE]. 

A stronger justification comes, however, from its effectivness in reproducing 
the numerical results, as detailed in |24) . From the simulations it is seen that 
the relaxation time of the local field tu = M'/T', rescaled by the mean collision 
time, increases with the packing fraction and with the inelasticity, as expected. 
At high densities it appears that r' ~ 1/0, and T tr ~ T g ~ Tg , likely due to 
stronger correlations among particles. At large <f> we observe T tr > Tj~, con- 
sistent with a smaller dissipation for correlated collisions. Model @ predicts 
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C = fc(t) and R = f R (t) with 

fc(R) = e' gt [cos(wt) + a C (R) sin(wt)] . (3) 

The variables <?, oj, ac and an are known algebraic functions of T E , T g , F', M' 
and T b . In particular, the ratio a c /a R = [T g - f2{T b - T g )]/[T g + Q(T b - T g )}, 
with fl = r E /{{r' + r E )(r E M'/M - r')). Hence, in the elastic (T g -> T b ) as 
well as in the dilute limit (J" — > oo), one gets ac — an and recovers the EFDR 
C(t) = R(t). Such predictions are all verified in numerical simulations [23] ■ In 
particular Fig. Q] depicts correlation and response functions in a dense case 
(elastic and inelastic): symbols correspond to numerical data and continuous 
lines the analytical curves. In the inelastic case, deviations from EFDR R(t) = 
C(t) are observed. In the inset of Fig. [1] the ratio R{t)/C(t) is also reported. 
It is important to notice that the main responsibility for the breakdown of the 
EFDR is the coupling between V and U, indeed Eq. ^ can be expressed in a 
different way: R(t) = aC(t) + b(V{t)U{0)) with a = [1 - (T g - T b )C2 a /r'] and 
b = {T g — T b )f2 b , where £2 a and f2 b are known functions of the parameters. At 
equilibrium or in the dilute limit the EFDR is recovered. 
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Fig. 1 (Color online). Left: correlation function C(t) (black circles) and response function 
R(t) (red squares) for a = 1 and a = 0.6, at <j> = 0.33. Continuous lines show curves obtained 
with Eqs. i3l . Inset: the ratio R(t)/C(t) is reported in the same cases. Right: Check of the 
fluctuation relation ([5} in the system with a = 0.6 and 4> = 0.33. Inset: collapse of the 
rescaled probability distributions of St at large times onto the large deviation function. 



An important independent assessment of model (J2J comes from the study 
of the fluctuating entropy production |25j which quantifies the deviation from 
detailed balance in a trajectory. Given the trajectory in the time interval 
[0,i], {V(s)Y , and its time-reversed {XV(s)}q = {— V(t - s)}q, the entropy 
production for our model takes the form ,20! 



This functional vanishes exactly in the elastic case, a = 1, where equipartition 
holds, T g — Tb, and is zero on average in the dilute limit, where (VU) = 0. 
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Formula (jU) reveals that the leading source of entropy production is the en- 
ergy transferred by the "force" PeU on the tracer, weighed by the difference 
between the inverse temperatures of the two "thermostats" . Therefore, to mea- 
sure entropy production, we need to measure the fluctuations of U: a possible 
choice is a local average of particles' velocities in a circle of radius I + R cen- 
tered on the tracer. Details on how to choose in a reliable way the correct I are 
given in [23]. Following such procedure, in the case = 0.33 and a — 0.6, we 
estimate for the correlation length I ~ 9r ~ 61q. Then, measuring the entropy 
production from Eq. Q along many trajectories of length t, we computed the 
probability P(St = x) and compared it to P(St = —%), in order to verify the 
Fluctuation Relation [Till 12 II 13] 

l0g P{E t = - x ) = X - (5) 

In the right frame of Fig. [T] the results of this comparison are reported. The 
main frame confirms that at large times the Fluctuation Relation ((SJ) is well 
verified within the statistical errors. The inset shows the collapse of log P(E t )/ t 
onto the large deviation rate function for large times. Notice that - in for- 
mula ([Jl - a wrong evaluation of the weighing factor (l/T g — 1/T&) or of the 
"energy injection rate" rEU{t)V{t) in Eq. ([J} could produce a completely 
different slope in Fig. [1] (right frame) . 

To conclude this paper, we stress that velocity correlations (V(t)U(t')) 
between the intruder and the surrounding velocity field are responsible for 
both the violations of the EFDR and the appearance of a non-zero entropy 
production, provided that the two fields are at different temperatures. We also 
mention that larger violations of EFDR can be observed using an intruder 
with a mass equal or similar to that of other particles |18] , with the important 
difference that in such simple "Langevin-like" model for the intruder's 

dynamics is not available. 
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